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ABSTRACT 

Stochastic acceleration of electrons and protons by waves propagating parallel 
to the large scale magnetic fields of magnetized plasmas is studied with emphasis 
on the feasibility of accelerating particles from a thermal background to rela- 
tivistic energies and with the aim of determining the relative acceleration of the 
two species in one source. In general, the stochastic acceleration by these waves 
results in two distinct components in the particle distributions, a quasi-thermal 
and a hard nonthermal, with the nonthermal one being more prominent in hot- 
ter plasmas and/or with higher level turbulence. This can explain many of the 
observed features of solar flares. Regarding the proton to electron ratio, we find 
that in a pure hydrogen plasma, the dominance of the wave-proton interaction 
by the resonant Alfven mode reduces the acceleration rate of protons in the in- 
termediate energy range significantly, while the electron-cyclotron and Whistler 
waves are very efficient in accelerating electrons from a few keV to MeV energies. 
The presence of such an acceleration barrier prohibits the proton acceleration 
under solar flare conditions. This difficulty is alleviated when we include the 
effects of 4 He in the dispersion relation and the damping of the turbulent waves 
by the thermal background plasma. The additional 4 He cyclotron branch of the 
turbulent plasma waves suppresses the proton acceleration barrier significantly 
and the steep turbulence spectrum in the dissipation range makes the nonthermal 
component have a near power law shape. The relative acceleration of protons 
and electrons is very sensitive to a plasma parameter a = u pe /fl e , where u pe and 
Q e are the electron plasma frequency and gyro-frequency, respectively. Protons 
are preferentially accelerated in weakly magnetized plasmas (large a). The for- 
malism developed here is applicable to the acceleration of other ion species and 
to other astrophysical systems. 
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1. INTRODUCTION 

One of the important questions in acceleration of cosmic particles is the fractions of 
energy that go into acceleration of electrons and protons (and other ions). In this paper, 
we investigate this question for acceleration by plasma wave turbulence, a second order 
Fermi acceleration process, which we call stochastic acceleration (SA). The theory of SA 
has received little attention in high energy astrophysics except in solar flares where it has 
achieved significant successes during the past few years. The turbulence or plasma waves 
required for this model are presumably generated during the magnetic reconnection which 
energizes the flares. The first application of SA was to the acceleration of protons and other 
ions to explain the observed nuclear gamma ray lines from solar flares (see e.g Ramaty 1979; 
Miller & Roberts 1995). Combining with the nuclear reaction rates (Ramaty, Kozlovsky & 
Lingenfelter 1975; 1979; see also Kozlovsky, Murphy & Ramaty 2002) and a magnetic loop 
model, Hua, Lingenfelter and Ramaty (1987a; 1987b; 1989) showed that the SA model can 
provide natural explanations for the many observed features in the 1 to 7 MeV range. Later 
this model was also investigated in the acceleration of electrons in several studies (Miller 
& Ramaty 1987; Bech, Schlickeiser & Steinacker 1989; Miller, LaRosa & Moore 1996; Park 
& Petrosian 1995 and 1996), and the first quantitative comparison of predictions of this 
model with the observed hard X-ray (10 to 200 keV) spectra in some solar flares was carried 
out by Hamilton & Petrosian (1992). With a more detailed modeling, Park, Petrosian & 
Schwartz (1997) showed that the SA of electrons by some generic turbulent plasma waves 
can reproduce the many spectral breaks observed over a broad energy range, from tens of 
keV to ~ 100 MeV, in the so-called electron dominated flares via the bremsstrahlung process 
(Rieger, Gan & Marschhauser 1998; Petrosian, McTiernan & Marschhauser 1994). 

The strongest evidence supporting the SA model comes from the YOHKOH discovery 
of impulsive hard X-ray emission from the top of a flaring loop, in addition to previously 
known emission from loop foot-points (FPs for short. Masuda et al. 1994; Masuda 1994). 
The presence of the loop-top (LT for short) emission requires temporary trapping of the 
accelerated electrons at the top of the loop where the reconnection is taking place. The 
turbulence required for SA will naturally accomplish this by repeated scatterings of the 
particles (Petrosian & Donaghy 1999). More importantly, an analysis of a larger sample of 
YOHKOH flares (Petrosian, Donaghy & McTiernan 2002) has shown that the LT emission 
is a common property of all flares, and a preliminary investigation of RHESSI data appears 
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to confirm this picture (Liu et al. 2003). Finally, a third and equally important evidence 
in support of the SA model comes from the spectra and relative abundances of the flare 
accelerated protons and other ions observed at 1 AU from the Sun (Mazur et al. 1992; Reames 
et al. 1994; Miller 2003). These several independent lines of arguments have established the 
SA as the leading model for solar flares. This may be also true in many other astrophysical 
nonthermal sources. Thus, a more detailed investigation of the SA model and its comparison 
with observations are now fully warranted. 

In particular, the SA of electrons on the one hand and protons and other ions on the 
other are investigated separately, a unified treatment and comparison with the total non- 
thermal radiative signatures of all species have not been carried out yet. The purpose of 
this investigation is to obtain the relative acceleration of electrons and protons from the 
thermal backgrounds of solar flare plasmas with the same spectrum of turbulence. We will 
present some general results of the model and qualitative comparisons with observations. 
More detailed comparisons with observations and the acceleration of other ions, such as the 
anomalous overabundance of the flare accelerated 3 He, will be addressed in subsequent pa- 
pers. Specifically, we will address the energy partition between the flare accelerated electrons 
and protons. Observationally, in some flares, or during the earlier impulsive phase of most 
flares, there is little evidence for gamma-ray lines and therefore proton acceleration. These 
are called electron dominated cases. In the majority of solar flares the energy partition 
favors electrons but there are a significant fraction of flares where more energy resides in 
protons than in electrons in their respective observed energy bands. The ratio of energy of 
the observed electrons (with > 20 keV range) to that of protons (with > 1 MeV range) in 
solar flares varies approximately from 0.03 to 100 (see e.g. a compilation by Miller et al. 
1997). In what follows, we will use solar flare plasma conditions but the formalism described 
here will be applicable to other astrophysical sources. 

In § 2 we describe the general theory of SA and argue that in most cases, the Fokker- 
Planck (F-P) equation can be reduced to the diffusion-convection equation with the corre- 
sponding coefficients given by pitch angle averaged combinations of the F-P coefficients. In 
§ 3, we study the resonant interactions of electrons and protons with parallel propagating 
waves in a pure hydrogen plasma, and calculate the resultant F-P coefficients and accelera- 
tion parameters for interactions with a power law turbulence spectrum of the wavenumber. 
The new and surprising result here is that the proton acceleration is suppressed by a barrier 
in its acceleration rate in the intermediate energy range. This is what is required by obser- 
vations of electron dominated cases, but as we will show this barrier is too strong and makes 
the acceleration of protons unacceptably inefficient relative to the electron acceleration, ex- 
cept for in very weakly magnetized plasmas. In § 4 we point out that this difficulty can 
be alleviated by a more complete description of the dispersion relation which includes the 
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effects of helium ions and by an inclusion of the effects of the thermal damping of the turbu- 
lence at high wavenumbers. The presence of an appropriate amount of fully ionized helium 
introduces an extra wave branch which lowers the barrier and the thermal damping steepens 
the turbulence spectrum toward high wavenumbers, making the acceleration of electrons and 
protons more in agreement with observations. The results presented here are summarized in 
§ 5 and their applications to solar flare observations are discussed qualitatively. Some useful 
approximate analytical expressions for the interaction rates are presented in the appendix. 



2. GENERAL THEORY OF STOCHASTIC ACCELERATION 

In this section, we present the general theory of SA and show that in most astrophysical 
situations, the diffusion-convection equation is adequate to address the particle acceleration 
processes. 



2.1. Fokker-Planck Equation 



The study of SA in a magnetized plasma starts from the collisionless Boltzmann-Vlasov 
equation and the Lorentz force (Schlickeiser 1989). In the quasilinear approximation, it can 
be treated by the F-P equation (e.g. Kennel & Engelmann 1966): 
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where the wave-particle interaction is parameterized by the F-P coefficients Dij[i,j e (fJ>,p)]- 
Here f(t, s,p, //) is the gyro-phase averaged particle distribution, s, v, /i and p are the spatial 
coordinate along the field lines, the velocity, the pitch angle cosine, and the momentum of 
the particle, respectively. The energy loss (minus systematic energy gains, if any) processes 
are accounted for by p\ and S is the source function. 

For weak turbulence (8B as is the case for solar flares, the F-P coefficients can 

be evaluated by assuming that the particles and waves are coupled via a resonant process. 
The acceleration of particles at a given energy is then dominated by interactions with certain 
specific wave modes, e.g. the Alfven or Whistler waves. For a study of acceleration in a 
narrow energy band it is usually sufficient to consider waves in a narrow frequency range 
(Miller & Ramaty 1987). In order to address the energy partition between electrons and 
ions, however, one has to calculate the particle acceleration over the whole energy range. For 
example, the Alfven waves can efficiently accelerate ions but not nonrelativistic electrons. 
Models with pure Alfvenic turbulence are not adequate to address the energy partition of 
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accelerated particles in solar flares and many other astrophysical plasmas, especially for 
the acceleration of particles from a thermal background. For the acceleration of such low 
energy particles interactions with turbulent plasma waves extending over a broad range of 
wavenumbers (and frequencies) must be considered. We consider here interactions with 
a broad spectrum of plasma waves propagating along a static background magnetic field. 
The interactions of parallel propagating waves with electrons are described fully by Dung 
& Petrosian (1994) (DP94) and Pryadko k Petrosian (1997) (PP97) (see also Steinacker & 
Miller 1992). We will use their formalism and evaluate the relative rates of interaction and 
acceleration of electrons and protons in cold but fully ionized plasmas. 



2.2. Dispersion Relation and Resonance Condition 

Waves propagating parallel or anti-parallel to the large scale magnetic field in a uniform 
plasma have two normal modes that are polarized circularly (Sturrock 1994). Because their 
electric fields are perpendicular to their corresponding wave vectors, the waves are also 
referred to as transverse waves. The dispersion relation for these waves in a cold plasma is 
(see DP94) 



{kef = uo 



(2) 



u{u + flj) 

where u P i = a/ 47m i g i 2 / 'm, and f2j = (qiB ) / (tojc) are respectively the plasma frequencies 
and the nonrelativistic gyro- frequencies of the background particles (with charges masses 
rrii and number densities rij). B stands for the large scale magnetic field, c is the speed of 
light, uj and k are the wave frequency and wavenumber, respectively. 

One of the important parameters characterizing a magnetized plasma is the ratio of the 
electron plasma frequency to the electron nonrelativistic gyro-frequency: 

a = uj pc /n c = 3.2(n e /10 10 cnr 3 ) 1/2 (A)/100 G)" 1 , (3) 

where fl e = (eB )/(m e c) and e and m c are the elemental charge unit and the electron mass, 
respectively. The value of a is small for a strongly magnetized plasma. 

A particle with a velocity f3c (Lorentz factor 7) and a pitch angle cosine /x interacts 
most strongly with waves satisfying the resonance condition: 

u - k\\Pn = — , (4) 

where n is for the harmonics of the gyro-frequency (not to be confused with the background 
particle number densities n,), uj and k\\ are the wave frequency and the parallel component 
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of the wave vector in units of ft c and ft e /c, respectively (we will use these units in the 
following discussion unless specified otherwise and in our case k\\ — k and n = —1), uJi = 
qifrie/emi is the particle gyro-frequencies in units of ft c ; u e = — 1 for electrons and u p = 
5 = m c /m p for protons, where m p is the proton mass (for more details see also DP94). One 
notes that low energy particles mostly resonate with waves with high wavenumbers and only 
relativistic particles interact with large scale waves with low frequencies. The resonant wave- 
particle interaction can transfer energy between the turbulence and particles with the details 
depending on the particle distribution and the spectrum and polarization of the turbulence. 



2.3. Fokker-Planck Coefficients 

The evaluation of the F-P coefficients requires a knowledge of the spectrum of the turbu- 
lence. Following previous studies (DP94; PP97) we will first assume a power law distribution 
of unpolarized turbulent plasma waves. For unpolarized turbulence, the amplitudes of the 
waves only depend on k. Then we have £{k) — (q — l)£ t ot^m7n^~ 9 ^ or ^ > ^ min (^ e - a l ar S e 
scale cutoff), where the turbulence spectral index q > 1. For a given turbulence energy 
density £ to t, k min , presumably larger than the inverse of the size of the acceleration region, 
determines the maximum energy that the accelerated particles can reach and the character- 
istic time scale of the interaction. The general features of this situation have been explored 
in the papers cited above. For the sake of completeness, we briefly summarize the key results 
here. 

The F-P coefficients can be written as 

. _ 2 . n ( /i/i(l — Xj) 2 , for ab = fifx; 

Dab = — — — X( k i) { W x A l ~ x j)> for ab = W\ (5) 
Tpi ^ j=i [ p 2 x 2 , for ab = pp, 



where 



The sum over j is for the resonant interactions discussed in the previous section. The 
characteristic interaction time scale for each of the charged particle species is r pi = t p /uj 
with that for electrons given by (see DP94): 

£tot 
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T- 1 ~ -ft 
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(9 ~ 1)A££ • (7) 



In general, the F-P coefficients have complicated dependence on the turbulence spectral 
index q, the plasma parameter a and the energy and pitch angle of the particles. The exact 
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solution of the full F-P equation is a difficult task. Fortunately under certain conditions 
considerable simplifications are possible. These conditions are defined by the relative values 
of the three F-P coefficients. The pitch angle change rate of the particles is proportional to 
D^, while the momentum or energy change rate is proportional to D pp /p 2 . As evident from 
equation (5) the behavior of D^ p /p is intermediate between the two. 



2.4. Diffusion-Convection Equation 

The relative values of the F-P coefficients determine the type of approximations that 
can be used for solving the F-P equation. We now show that for most conditions reasonable 
approximations lead to the well known transport equation (eq. [10]). In order to justify 
these approximations it is convenient to define two ratios of the coefficients: 

R^p) = -g^, (8) 

We will show in the following sections, for most energies and pitch angles both Ri and |i?2| ^ 
1, which means that ^> D pp /p 2 . Under these conditions the particles are scattered 
frequently before being significantly accelerated and the accelerated particle distribution is 
nearly isotropic. Then the pitch angle averaged particle distribution function F(s,t,p) = 
0.5 J*_ 1 d///(/i, s, t,p) satisfies the well known diffusion-convection equation (see e.g. Kirk, 
Schneider & Schlickeiser 1988; DP94; PP97). 

In this study we are interested in the relative acceleration of electrons and protons which 
is not sensitive to the detailed geometry or the inhomogeneities of the source. Therefore we 
can assume a homogeneous and finite (size L) source, or alternatively confine our discussion 
to spatially integrated spectra. In this case we can treat the spatial diffusion or advection 
of the particles by an energy dependent escape term. Then the above mentioned equation 
is reduced to 

dN d 2 d N 

;(D EE N) + ^-[(E h -A)N]- — + Q, (10) 



dt dE 2X arj ' <9£ LV ' J T r 



CSC 



where E — (7 — l)mjC 2 is the particle kinetic energy, N(t,E)dE = Anp 2 dp ft F(s,t,p)ds, 



El describes the net systematic energy loss, and Q(t, E) = 0.5 d/i J Q L S(s, /J,, t, E)ds is the 
total injection flux of particles into the acceleration region. Dee describing the diffusion in 
energy is related to D pp and defines the acceleration time, and T esc is related to the scattering 
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time r sc : 

T esc = (L 2 /v 2 )/r sc , r sc = \ f d^ (1 ~^ )2 « L/u , (11) 
r ac = S 2 /£) EB , D EE = — J ^ dfiD^Ri - R\) . (12) 

Note that equation (10) describes the energy diffusion with two terms, D EE and the direct 
acceleration rate: 

1 d(3 7 2 D EE dD EE D EE 2- 1 ~ 2 

A{E) = W 2 ^e~ = ^e~ + ITT+r 1 ■ (13) 

There are several important features in the diffusion coefficients which we emphasize 

here. 

1) The first is that in the extremely relativistic limit the diffusion coefficients (and their 
ratios) for protons and electrons are identical and assume asymptotic values such that both 
of the ratios are much less than one. Therefore equations (10), (11) and (12) are valid. 
(Strictly speaking, this is not true for very strongly magnetized plasmas a < 5 l l 2 where one 
gets R\ ~ \R 2 \ ~ 1. See eq. [5]) 

2) The second is that at low energies, as pointed out by PP97, R\ and R\ are not 
necessary less than one, especially for plasmas with low values of a. In the extreme case of 
Ri 3> I-R2I 3> 1, three of the four diffusion terms in equation (1) can be ignored. Again, if 
we assume a finite homogeneous region, or integrate over a finite inhomogeneous source, the 
resultant equation becomes similar to equation (10). Now because of the lower rate of pitch 
angle scatterings, the escape time may be equal to the transit time T esc ~ L/(vfi), the other 
transport coefficients D EE and (and consequently the accelerated particle spectra) may 
depend on the pitch angle, and the assumption of isotropy may not be valid. However, as can 
be seen in the next section (Figures 5 and 6), these coefficients change slowly with //, except 
for some negligibly small ranges of //, so that the expected anisotropy is small. In addition, at 
lower energies Coulomb scatterings become increasingly important and can make the particle 
distribution isotropic. In many cases, especially for plasmas not completely dominated by the 
magnetic field (i.e. for a > 1) one can neglect the small expected anisotropy and integrate 
the equation over the pitch angle, in which case the transport equation becomes identical to 
equation (10) except now 

T esc = L/V2K<r sc ~<l/^>, (14) 
2p 2 

Tac = r 1 a n rv (15) 

J_i d^D pp {fi) 
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where "<>" denotes averaging over the pitch angle. 

3) It is easy to see that one can combine the above two sets of expressions (eq. [11] -[15]) 
for the acceleration rates (or time scales) and the escape times at the nonrelativistic and 
extremely relativistic cases as 



The first expression in equation (17) is valid at low values of E and a and the second at 
higher energies and in weakly magnetized plasmas. However, it turns out that at extremely 
relativistic energies and in weakly magnetized plasmas (a > 1), independent of other con- 
ditions, R\ <C R\ and the first expression can be used. These expressions and equation 
(10) then describe the problem adequately for most purposes in high energy astrophysics, in 
particular for solar flares, the focus of this paper. 

4) Finally, in certain cases, especially in the intermediate energy range the quantity 
Ri — R\ appearing in equations (12) and (17) can be small. The acceleration rate can be 
reduced dramatically when both R\ and |i? 2 | are much less than one and R\ ~ R\. From 
the definitions of these ratios and expressions for the F-P coefficients (eqs. [8], [9] and [5]) it 
is clear that if there were only one resonant interaction one would have R\ = R\ and there 
would be no acceleration. Thus, strictly speaking the use of equation (10) with interactions 
involving only one wave mode (say the Alfven) is incorrect. However, as we will show in § 3.1 
there are always at least two resonant interactions in unpolarized turbulence, in which case 
Ri 7^ R\ so that the acceleration rate is finite. But if one of the interactions is much stronger 
than the others, R\ — R\ can be small. In the next section, we will show some examples 
where this is true (Figure 5) and that this happens at the intermediate values of energy 
(Figure 6). The acceleration rate is then reduced greatly. The much lower acceleration 
rate at the intermediate energies compared to the higher rates in the nonrelativistic and 
extremely relativistic limits introduces an acceleration barrier. As we shall see that in the 
intermediate energy range the behaviors of protons and electrons are quite different and a 
much stronger acceleration barrier appears for protons. 




(16) 



and 




(17) 
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2.5. Loss Rate 

To determine the distributions of the accelerated protons and electrons by solving equa- 
tion (10) with the above formalism, in addition to the transport coefficients Dee, A and 
T esc , we need to specify the loss term E L . For electrons the loss processes are dominated by 
Coulomb collisions at low energies and by synchrotron losses at high energies: 

E Le = 4r 2 m e c 3 [7m c In A//3 + S 2 /3 2 7 2 /9m e c 2 ] , (18) 

where r = 2.8 x 10~ 13 cm is the classical electron radius and In A = 20 is a reasonable 
value in our case (See Leach 1984). The ion losses in a fully ionized plasma are mainly due 
to Coulomb collisions with the background electrons and protons (Post 1956; Ginzburg & 
Syrovatskii 1964). For electron-ion collisions, we have 

for (5 < (3 Te ; 
for 1 > (3 > (3 Te ; 

( 19 ) 

for mi/m c ^> 7 ^> 1 ; 
for 7 > rrii/me , 

where f3^ e = {"ik^T c /m c c 2 ) 1 ' 2 is the mean thermal velocity of the background electrons in 
units of c and k\> is the Boltzmann constant. For proton-ion collisions, which are important 
for ions with even lower energies, we have (Spitzer 1956) 

E Li = 47rr 2 m e c 3 n p (g i /e) 2 (m c /m p )/3- 1 lnA for f3 > /3 Tp - (20) 

where (3t p = (3A;bT p /m p c 2 ) 1//2 is the mean thermal velocity of the background protons in 
units of c. These loss processes dominate at different energies and we can define a loss time 

Tloss = E/E L . 

2.6. Steady State Solution and Normalization 

We will use the impulsive phase conditions of solar flares for our demonstration. In this 
case, we can assume that the system is in a steady state because the relevant time scales 
are shorter than the dynamical time (the flare duration). We also assume the presence of 
a constant spectrum of turbulence. We are interested in the acceleration from a thermal 
background plasma, therefore, a thermal distribution is assumed for the source term Q. 
As described above, equation (10) may not be valid at low (keV) energies where R\ 3> I. 
However, for solar flare conditions and in the keV energy range, Coulomb scatterings become 



f 2 v /6/^ 2 /^ 3 lnA 
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important (D^° ul > D™ ve , see Hamilton & Petrosian 1992). In this case Ri < 1 and the 
particle distribution will be nearly isotropic at all energies. We therefore calculate the 
acceleration rate with the second expression of equation (17) and solve equation (10) to get 
the distributions of the accelerated particles over all energies. 

To appreciate the relevant physical processes, one can compare the acceleration time 
with the escape and the loss time. We are mostly interested in the energy range above 
the energy of the injected particles. So the source term is not as important in shaping 
the spectrum as the other terms. In the energy band where the escape and loss terms 
are negligible, from the flux conservation in the energy space, one can show that AN — 
&(DeeN) /&E = constant. On the other hand, when the acceleration terms are negligible, 
no acceleration occurs. When the escape time becomes much shorter than the acceleration 
time and both of them are much shorter than the loss time, particles escape before being 
accelerated. This results in a sharp cutoff in the particle distribution at the energy where 
T esc ~ E/A(E) ~ E 2 /Dee- When the escape time is long and the loss time is much 
shorter than the acceleration time, one would then expect a quasi-thermal distribution for 
the Coulomb collisional losses (Hamilton & Petrosian 1992) and a sharp high energy cutoff 
for the synchrotron losses (Park, Petrosian & Schwartz 1997). Power-law distributions can 
be produced only in energy ranges where the loss term is small and the acceleration and 
escape times have similar energy dependence. 

The normalization of the steady state particle distributions is determined by their rates 
of acceleration, escape and injection. The injection rates depend on the geometries of the 
reconnection and the turbulent acceleration site and on possible contributions of the charged 
particles to reverse currents which must exist when a net charge current leaves the acceler- 
ation site. A more detailed time dependent treatment is required to determine the relative 
normalization. This is beyond the scope of the paper and will be dealt with in the future. 
Here we concentrate on the relative shapes of the electron and proton spectra in the LT and 
FP sources. We will assume that the injection flux J QdE = 1 s^ 1 cm -2 for both electrons 
and protons (see also § 5). In the steady state this is equal to the flux of the escaping particles 
TV*"* = J °° Nlt^E) /T esc (E)dE . Since the escaping particles lose most of their energy at the 
FPs, instead of N esc (E) = N LT (E)/T esc (E) we will show the effective particle distribution 
for a thick target (complete cooling) FP source, which is related to the corresponding LT 
distribution 7V L t via (Petrosian & Donaghy 1999): 




(21) 



3. APPLICATION IN COLD HYDROGEN PLASMAS 



In this section we describe the relative acceleration of electrons and protons in cold, 
fully ionized, pure hydrogen plasmas. This is an approximation because all astrophysical 
plasmas contain some helium and traces of heavy elements. Ignoring the effects of helium 
(trace elements will, in general, have no influence on the following discussion) and adopting 
a turbulence spectrum of a single power law of the wavenumber simplify the mathematics 
and allow us to demonstrate the differences between the acceleration rates of electrons and 
protons more clearly. Moreover, in some low temperature plasmas, most of the helium may 
be neutral and not be involved in the SA processes. The results presented here will be a 
good approximation. Pure hydrogen plasmas can also be realized in terrestrial experiments 
to test the theory. The formalism can also be easily generalized to the case of electron- 
positron plasmas and to more complicated situations. In the next section we will present 
our results for plasmas including about 8% by number of helium and for turbulence with a 
more realistic spectrum. 



3.1. Dispersion Relation and Resonant Interactions 

In a pure hydrogen plasma, equation (2) reduces to (PP97) 

^ = 1 - a2{1 + 6) (22) 
oo 2 (w-l)(w + <5)' 1 1 

and the Alfven velocity in units of c is given by (3\ = 5 1 / 2 / a. (For e ± pair dominated plasmas 
5 = 1). 

Figure la depicts the normal modes of these waves, which compose four distinct branches. 
From top to bottom, we have the electromagnetic wave branch (EM; long dashed), electron- 
cyclotron branch (EC; dot-dashed), proton-cyclotron branch (PC; dotted), and a second 
electromagnetic wave branch (EM'; short dashed), respectively. The lower panel is an en- 
largement of the region near the origin. The positive and negative frequencies mean that the 
waves are right- and left-handed polarized, respectively, where the polarization is defined 
relative to the large scale magnetic field (Schlickeiser 2002). Figure lb depicts the group 
velocities f3 g = duj/dk of these waves. One may note that the signs of the phase velocity 
(3 ph = u/k and the group velocity of a specific wave mode are always the same. 

In Figure la, the two solid straight lines depict equation (4) for an electron (upper) and 
a proton (lower) with f3 = 0.5 and fi = 0.25. The intersections of these lines with the wave 
branches satisfy the resonance condition. The electron interacts resonantly at the indicated 
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point with the EC branch and at another point with the PC branch at a high negative 
wavenumber that lies outside the figure. The proton, on the other hand, not only resonates 
with one PC wave, but also with three EC waves (only two of which are seen in the lower 
panel of the figure). As we shall show below, the fact that certain protons can resonate 
with more than one EC wave has significant implications for the overall proton acceleration 
process. 
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Fig. 1. — a) Left panel: Dispersion relation of parallel propagating waves in a cold pure hydrogen 
plasma with a = 0.5. The bottom panel is an enlargement of the region around the origin. The 
curves, from top to bottom, describe the EM (long-dashed), EC (dot-dashed), PC (dotted) and EM' 
(short-dashed) waves. The upper and lower solid lines give respectively the resonance conditions 
for electrons and protons with v = 0.5c {(3 = v/c) and fi = 0.25. Resonant interactions occur at the 
points where these lines cross the curves which depict the waves, b) Right panel: Same as the 
left panel but for the group velocity /3 g = dui/dk versus the wavenumber k. The line type remains 
the same for each wave branch. Negative group velocities mean that the energy fluxes of the waves 
are in the direction anti-parallel to the large scale magnetic field. 
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3.2. Critical Angles and Critical Velocities 

In general one expects four resonant points. However, for a given particle velocity or 
energy, at critical angles, where the group velocities of the waves are equal to the parallel 
component of the particle velocity, the number of resonant points can change from four to 
two or vice verse. Figure 2 shows the velocity dependence of the critical angles for electrons 
(upper panels) and protons (lower panels) in plasmas with a = 0.5 (left panels) and a = 0.1 
(right panels). (The results for electrons are the same as those given by PP97.) Both particles 
have at least two resonant interactions (one with the PC and one with the EC branch except 
for /i — where electrons interact with two EC waves and protons interact with two PC 
waves) . 

Electrons with a large /j, can have two additional resonances with the EM branch and 
those with a small \x have two additional resonances with the EC branch. The two regions 
with four resonances grow with decreasing a and shrink as a increases. For larger values of 
a the interaction is weaker because for large ranges of velocities and pitch angles electrons 
interact with only two waves (e.g. the interactions with the EM branch disappear for a > 1. 
See Figures 2 and 3 ). But as a approaches zero the region with two wave interactions 
diminishes and the two curves for the critical angles merge into one, satisfying the relation 

/i cr = ( 7 -l)//3 7 , for a^O. (23) 

In the case there are always four resonances and the total interaction is strong at all energies. 

Protons have a similar, but slightly more complicated, behavior. As fi increases one 
obtains interactions with 1EC+3PC, 1EC+1PC, 3EC+1PC and back to 1EC+1PC waves. 
With the decrease of a, the upper two regions diminish, while the lower portions increase in 
size. Protons can also be accelerated by the EM' waves but this only occurs in more highly 
magnetized plasmas (a < 5 1 ^ 2 /2 ~ 0.012) as compared with the interactions of electrons with 
the EM branch. At such low values of a a region with four interactions (1EC+1PC+2EM') 
appears in the upper portion of the /i — (5 plane and its lower boundary eventually merges 
with the lower curve for /i cr as a approaches zero. Just like electrons the critical angle is 
given by equation (23). In this limit, particles are basically exchanging energy with the 
Poynting fluxes of the electromagnetic waves. 
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Fig. 2. — Left Panel: Velocity dependence of the critical angles in a plasma with a = 0.5 for 
electrons (top) and protons (bottom). At small pitch angles, i.e. fj, — ► 1, electrons can resonate with 
high frequency electromagnetic waves of the EM branch (region labeled with"2EM+lEC+lPC"), 
while energetic protons mainly interact with the Whistler and Alfven waves (region labeled 
with"3EC+lPC"). Right Panel: Same as above but for a = 0.1 where interactions of protons 
with the Whistler waves start at a higher energy. See text for details. 
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These behaviors can also be seen in Figure 3, where instead of /x cr we plot what one 
may call the critical velocities as a function of a for two values of fi. Note that in the proton 
panel with \x = 1.0 there is a small region with a < 0.012 where there are four resonances 
including two with the forward-moving left-handed polarized electromagnetic waves from the 
EM' branch. Protons will not resonate with the electromagnetic waves for larger values of 
a. In general, we have similar patterns of transition between different regions caused by the 
electromagnetic branches in the /i — (5 space, except that the transitions for protons occur at 
a value of a which is lower than that for electrons by a factor of 5 1 / 2 . The main difference in 
the behaviors of electrons and protons resides in their four resonant interactions with the PC 
and EC branches. Protons have two such regions where they resonate with "1EC+3PC" or 
"3EC+1PC", while electrons only have one with "3EC+1PC"; electrons never interact with 
more than one PC wave. This is where the above scaling symmetry of a between protons 
and electrons is broken. 
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Fig. 3. — Same as Figure 2 but depicting the dependence of the critical velocities on the plasma 
parameter a. Combining this with Figure 2, one can tell the wave branches responsible for the 
critical velocities. 
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Low Energy Approximations: Because the acceleration of particles at low energies is of 
particular interest we present here some approximate analytic relations, which are derived 
in the appendix. 

The first is for the proton critical velocity curve dividing the region with two and four 
resonances (i.e. the middle curve in the lower left panel of Figure 3). As we will see in 
the following sections, at a given a the acceleration rate (eq. [17]) increases dramatically 
once protons attain the critical velocity or energy and enter the region with four resonant 
interactions. The pitch angle averaged acceleration rate also increases sharply above this 
energy. It will be useful to have a formula to estimate this critical velocity. We find the 
following approximate expression for this transition 



which is shown by the dotted curve in the low left panel of Figure 3 and agrees within 0.2% 
with the exact result for a > 0.05. 

The second approximation is for the critical angles of protons below the critical energy 
(velocity) and low energy electrons, most of which interact only with two waves with one 
dominating over the other. When this happens, the acceleration rate for the particles can 
be very small (see § 2.3). Only particles with very large pitch angles (fx ~ 0) have four 
resonances and significant contributions to the pitch angle averaged acceleration rate. The 
regions for this lie in the small areas below the lowest curves in Figure 2, which are barely 
visible for the proton and a = 0.5 case. As shown in the appendix, using the approximations 
of equations (Al) and (A2) for the dispersion relations, we can derive analytic expressions 
for the critical angles, which in the nonrelativistic limit give /x cr oc f3 2 oc E/mc 2 . Empirically, 
we find the following simple approximate expressions, as shown in Figure 4, agree with the 
exact results to better than ~ 10% in the indicated energy ranges: 



cr 



0.0012 + a 2 



0.06a 



or 




(24) 




5 1 l 2 E/m v c 2 , for protons; E < E CI , 
E/m c c 2 , for electrons; E < 60 keV. 



(25) 
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Fig. 4. — Comparison of the analytical expressions (eq. [25] solid curves) with the exact values 
of /i cr (dotted and dashed curves; in the low energy range) due to resonant interactions with the 
EC (left panel for electrons) or PC (right panel for protons. Note that the region to the left of 
EpOt 2 = 40 keV has an expanded scale) branch. The thin dashed (barely visible near the left axis 
and for a = 1.0) and dotted (for a = 0.1) curves in the right panel give the critical angles for waves 
obeying the Alfvenic dispersion relation u = —\k\flA-, which clearly give incorrect descriptions for 
the acceleration of low energy protons. Low energy electrons do not interact with the Alfven waves. 
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Here it should be emphasized that the commonly used approximation of accelerating 
protons by the Alfven waves with the dispersion relation u = —\k\^A f° r M < which 
is valid at relativistic energies (Barbosa 1979; Schlickeiser 1989; Miller & Roberts 1995), 
is invalid at low energies. This can be seen from the lower left panel of Figure 1, which 
shows clearly that the dispersion relation of the waves in resonance with low energy protons 
deviates far from the simple Alfvenic form. For the simple Alfvenic dispersion relation, 
nonrelativistic protons resonate with both a forward and a backward moving wave only if 
\n\ < fi cr = Pa(i — l)//^7- (Particles only interacting with one wave can not be accelerated; 
§ 2.4.) This critical angle is indicated by the thin dashed (barely visible near the vertical 
axis) and dotted curves in the right panel of Figure 4 for a = 1.0 and 0.1 respectively, which 
clearly overestimate by several orders of magnitudes the fractions of low energy protons that 
can be accelerated. The inefficiency of proton acceleration at intermediate energies combined 
with the increase of the interaction rate above E cr gives rise to the acceleration barrier to be 
described in §3.4. 

In most of the particle acceleration models, an injection process of high energy particles 
is postulated as an input. If the injected particles have an energy above E CT , it may be 
appropriate to use the Alfvenic dispersion relation to describe the waves. If the energy 
of the injected particles is low, as is the case under studied here, one must use the exact 
dispersion relation to calculate the acceleration of low energy particles. Although most of the 
turbulence energy is carried by waves with low wavenumbers, the acceleration of low energy 
particles is determined by waves with high wavenumbers (see eq. [4]), which can constrain 
the overall acceleration efficiency As discussed above, the Alfvenic dispersion relation for 
the waves will overestimate the acceleration efficiency of low energy protons significantly. 



3.3. Fokker-Planck Coefficients 

For a given power law spectrum of the turbulence, it is straightforward to calculate the 
F-P coefficients with equation (5). Figure 5 shows variation of these coefficients with /i at 
some representative energies and Figure 6 shows the variation of their inverses (i.e. times) 
with energy at different values of //, here a = 0.5 and q = 1.6. The discontinuous jumps 
occur at the critical values of /z cr and f3 cr described in § 3.2. The variations of the two ratios 
defined in equations (8) and (9) with energy at different pitch angles are shown in Figure 7. 
These results justify the discussion in § 2.4. 
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Fig. 5. — Pitch angle dependence of the F-P coefficients (absolute values) for protons of two 
different energies (8MeV: left panel and 1 GeV: middle panel), and 200 keV electrons (right panel). 
A, B, C, and D stand for (D pp — D ppi / D^)/p 2 , D pp /p 2 , \D p ^\/p, and D^, respectively. The plasma 
parameter a = 0.5 and the turbulence spectral index q = 1.6. Note that different coefficients are 
scaled differently and that for illustration, in the middle panel the region to the left of \i = 0.2 is 
expanded. 
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Fig. 6. — Energy dependence of the scattering, 1/D = D~^, and acceleration times, 1/B = p 2 /D pp 
or 1/A = p 2 / (Dp P — D 2 p /D^), in units of r p at different pitch angles for electrons (left panel) and 
protons (right panel). The plasma parameter a = 0.5 and the turbulence spectral index q = 1.6. 
Note that A = B for /i = 0. At very low energies > p 2 /D pp and the acceleration time 

T ac = P 2 /D pp . At intermediate and high energies r ac = p 2 /(D pp - D 2 /D p ^). 
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Fig. 7. — The ratios R\ (left panel) and R\ (right panel) as functions of energy for several different 
pitch angles for electrons (upper panels) and protons (lower panels). In general, R\ and |i?2| exceed 
unit at low energies and [i (and small values of a, see PP97). For protons, at certain ranges of E 
and /j,, R\ ss R\ <C 1, and the acceleration rate, D^{Ri — R\) can be very small, (see eq. [12]). 
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3.4. Barrier in the Proton Acceleration 

In the previous section, we showed that the pitch angle averaged acceleration rate is one 
of the dominating factors in the particle acceleration processes. The relative acceleration of 
protons and electrons therefore depends on the contrast of their acceleration times. Figure 
2 shows that at low energies particles with \i > /x cr only resonate with one PC and one EC 
wave: the EC wave dominates the PC wave for electrons while the reverse is true for protons. 
These particles have significant contributions to the pitch angle averaged acceleration rate 
for /i cr <C 1 (eq. [25]). Because the difference between the wavenumbers of the two waves 
interacting with protons is much larger than that for electrons, the resonant interaction is 
more strongly dominated by one of the resonant waves for protons than it is for electrons. The 
factor 1 — R\jR\ for protons can therefore be several orders of magnitude smaller than that 
for electrons at a given energy and pitch angle (Figure 6). Consequently, in the intermediate 
energies where Ri ~ R\ ~ 1, the pitch angle averaged acceleration time for protons has 
a more prominent increase than that for electrons. At still higher energies, particles with 
four resonances dominate the acceleration rate because R\ <C R\ for the interactions. For 
both electrons and protons, the new resonant wave modes come from the EC branch. In the 
relativistic limit, the acceleration is dominated by resonances with the Alfven waves and the 
interaction rates for electrons and protons become comparable. 

Figure 8 shows the pitch angle averaged acceleration (thick curves) and scattering (thin 
dashed curves) times in units of r p for electrons (lower curves) and protons (upper curves) in 
plasmas with a = 0.5 (left panel) and a — 0.1 (right panel). The turbulence spectral index 
q = 1.6 here. The acceleration times for both cases in equation (17) are plotted with the 
invalid segments plotted as dotted curves. We see that the pitch angle averaged acceleration 
times are much shorter than the corresponding scattering times for keV particles. The 
particle distributions can be anisotropic at these energies unless there are other scattering 
processes (e.g. Coulomb collisions). In the high energy range, the scattering time is always 
shorter than the corresponding acceleration time when f3^ < 1. The transitions where 
Ri ~ I-R2I ~ 1 (as indicated by the circles in the figure) occur between 10 2 to 10 3 keV, 
increase with the decrease of a and depend on the turbulence spectral index q as well. 
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Fig. 8. — Pitch angle averaged acceleration (thick curves) and scattering (thin dashed curves) 
times in units of r p in plasmas with a = 0.5 (a: left panel) and a = 0.1 (b: right panel). Here 
the turbulence spectral index q = 1.6. The acceleration times defined for both cases of equation 
(17) are plotted with the corresponding invalid segments indicated by the dotted curves. The upper 
two curves are for protons and the lower two are for electrons. The circles indicate the points of 
transitions between low and high energies where i?i ~ |i?2| ~ 1- The transitions of the electron 
acceleration times are quite smooth. The thin solid lines show schematically the transitions of the 
acceleration time for protons. The acceleration barrier (as indicated by the hatched areas) in the 
proton acceleration times are prominent. The sharp drop of the electron acceleration time with the 
decrease of energy for a = 0.1 is due to interactions with the EM branch (PP97). 
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There is clearly an acceleration barrier (as indicated by the shaded area) in the proton 
acceleration time. The thin solid line shows schematically the acceleration time of protons 
in the transition region. The sharp increase of the proton acceleration time at lower energies 
is caused by their low acceleration efficiency when the scattering rate already overtakes the 
acceleration rate as discussed above. The sharp drop of the proton acceleration time at a 
higher energy is due to interactions with the Whistler waves. Because protons with small 
pitch angles (// ~ 1) interact with the Whistler waves at the lowest energy and the interaction 
is very efficient, this energy corresponds to the critical energy E cr identified in equation (24). 

These characteristics are not true for electrons. In section § 3.2, we have shown that a 
small fraction of particles with ji < ji cv can resonate with four waves and |i? 2 | «C R\ for the 
interactions (Figure 5). Compared with protons, more electrons can be accelerated this way 
(eq. [25]). Because the acceleration of electrons with two resonances is very inefficient, the 
acceleration of this small fraction of electrons already dominates the electron acceleration 
processes where the scattering rate becomes comparable with the acceleration rate. At even 
higher energies, there is no extra wave modes which can enhance the electron acceleration 
processes. Consequently, electrons have a smooth acceleration time profile. 

Expressions for estimating the difference between electron and proton acceleration time 
are derived in the appendix. Briefly, because particles with [i < fi cr have significant contribu- 
tions to the acceleration in the low energy region (f3 <C 1), one can estimate the pitch angle 
averaged acceleration times with the approximate expressions for the critical pitch angles 
(eq. [25]): 



— = 2 



W +2 {AY~ q)/2 \]-^ f releC ; r ° nS; (26) 
\m e c 2 J yd 1 for protons, 



which is consistent with the numerical result within a factor of two. In the relativistic limit 
(7 ^> 1), particles interact with the Alfven waves and we find 

r ac q(q + 2)a 2 ( eV~ q ( E ^ 2 ~ q 



t p 45 \qij \m e c 



-,2 



(27) 



The difference between these two time scales at the critical energy (eq. [24]) gives an estimate 
of the height of the acceleration barrier. 



3.5. Application to Solar Flare Conditions 



Figure 9 shows a model of electron (thick curves) and proton (thin curves) acceleration 
in a strongly magnetized plasma. The LT size L = 10 9 cm, the temperatures of the injected 
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electrons and protons are the same k^T = 1 keV. The magnetic field and gas density are 400 
G and 4. Ox 10 9 cm -3 , respectively, i.e. a = 0.5 (see equation [3]). The relevant time scales are 
shown in the left panel, where we have defined the direct acceleration time r a = E/A, which 
is related to r ac (eqs. [12] and [13]). The corresponding accelerated particle distributions 
are shown in the right panel. The dotted and dashed curves show the LT and FP spectra, 
respectively. 

We note that for the above conditions the electrons can be accelerated to a few hun- 
dreds of keV while the proton acceleration is suppressed due to the acceleration barrier. The 
electron distribution steepens with the increase of energy because the escape time becomes 
shorter than the acceleration time (T esc < r a ). At low energies where Coulomb collisions dom- 
inate, the LT electrons have a quasi-thermal distribution. The solid curve gives the thermal 
distribution of the injected particles with arbitrary normalization. Due to the absence of 
acceleration at low energies, the steady state proton distributions are almost identical with 
the injected proton distribution. 

To produce a near power law electron distribution, as suggested by solar flare observa- 
tions, the escape time must be comparable with the acceleration time in the relevant energy 
band. Because the escape time of the nonrelativistic electrons always decreases with the 
increase of energy, we adopt a turbulence spectral index of q = 3 in the model. The plasma 
time t p = 1 s. Then we have the ratio of the turbulent wave energy density to the magnetic 
field energy density 8n£ tot /B^ = 4.4 x 10~ 10 A: m f n . To ensure that this ratio is much less 
than one so that the quasilinear approximation for the F-P treatment stays valid, one needs 
&min > 10~ 5 or an injection length of the turbulent waves of less than 10 6 cm (note that /c m j n 
is in units of Vl c /c ~ 0.24 cm -1 here). Otherwise, the turbulence spectrum must flatten at 
low k so that there is less energy content in long wavelength waves. 
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Fig. 9. — Left panel: Time scales for protons (thin curves) and electrons (thick curves) in a 
strongly magnetized plasma (a = 0.5) with a steep turbulence spectrum (q = 3). The direct 
acceleration times r a = E/A, which are related to the diffusion acceleration times r ac , are shown by 
the dotted curves. The solid curves are for t\ oss and the dashed curves are for T esc . Right panel: 
The corresponding distributions N * E 2 of the accelerated electrons (thick curves) and protons 
(thin curves). Here the dotted curves give the thin target LT particle distributions, the dashed 
curves indicate the effective thick target particle distributions at the FPs. The solid curve gives 
the injected thermal particle distribution with arbitrary normalization. One can see that because 
of the presence of the strong acceleration barrier (exceeding the range of the graph) protons are 
basically not accelerated. 



-30- 



In § 3.2 and § 3.4, we showed that the proton acceleration barrier moves toward lower 
energies with the increase of a. So in very weakly magnetized plasmas, this barrier can 
be close to the thermal energy of the injected particles and thus has little effect on the 
acceleration of protons. Protons can be accelerated efficiently in the case because their loss 
time is long. Figure 10 shows such a model, where B = 100 G and n = 10 11 cm" 3 . The size 
of the LT and the injected particle temperatures remain the same as those in the previous 
model. Because the turbulence spectrum is flat (q — 2), we have a pretty hard accelerated 
proton distribution below 1 MeV. Above this energy, there is a cutoff due to the dominance 
of the escape term over the acceleration terms. The accelerated electron distribution has a 
cutoff at less than 100 keV which is also due to the quick escape of electrons with higher 
energies from the acceleration site. At a few keV, both electron and proton distributions are 
quasi-thermal because of the dominance of Coulomb collisions. 

The above results show that electrons can be accelerated to very high energies by parallel 
propagating turbulent waves in pure hydrogen plasmas, but the presence of the acceleration 
barrier in the intermediate energy range makes the acceleration of protons very inefficient. 
Only in very weakly magnetized plasmas where the barrier is close to the background particle 
energy does the acceleration of protons become efficient. The required value of the plasma 
parameter a is above 10 which is much larger than that believed to be the case for solar flares. 
However, most astrophysical plasmas including solar flares are not made of pure hydrogen. 
They contain significant numbers of 4 He. These particles modify the dispersion relation used 
above. Abundances of elements heavier than He are too small to have a significant effect 
but 4 He with an abundance (by number) of about 8% can have important effects. 

To produce a near power law distribution of the accelerated particles, the index of 
the turbulence spectrum must be larger than 2 at high wavenumbers. If the turbulence is 
generated on a scale comparable to the size of the flaring loops, such a steep turbulence 
spectrum must flatten at low wavenumbers to ensure that the turbulence energy density is 
less than the energy density of the magnetic field, which is presumably the dominant source 
of energy for solar flares. Such a steepening of the turbulence spectrum at high k is expected 
if one includes the thermal damping effects of the waves with high wavenumbers. We will 
incorporate these effects in the following discussion. 
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Fig. 10. — Same as Figure 9 but for a model in a weakly magnetized plasma (a = 10). The model 
parameters are shown in the figure. Here the proton acceleration becomes more efficient than the 
electron acceleration because its barrier is close to the energy of the injected particles and does not 
affect the acceleration processes. 
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4. ACCELERATION IN HYDROGEN AND HELIUM PLASMAS 

We now repeat the derivation of the previous section for more realistic plasmas contain- 
ing e, p and a-particles. We assume a fully ionized H and 4 He plasma with the following 
relative abundances: electron/proton/o;-particle= 1/0.84/0.08. 



4.1. Dispersion Relation and Resonant Interactions 

The dispersion relation for such a plasma can be written as: 



(J 1 UJ 



1 , (1 - 2Y ne )S , Y He S 



+ 



u-l uj + 5 uj + 5/2 



(28) 



where the He abundance Y~He = 0.08. The other symbols are the same as those defined in § 
2.2. 

The inclusion of 4 He splits the PC branch into two: one covers the frequency range of 
to —5/2 and the other covers the frequency range of — upc to —5, where ujpc — (0.5 + Yn e )5 
is the lowest frequency of the branch (remember that the minus sign only indicates that 
the waves are left-handed polarized). We refer the former as a 4 He cyclotron branch (HeC) 
and the later as a modified proton cyclotron branch (PC) because at high wavenumbers 
they approach to the a-particle and proton cyclotron frequencies, respectively. Figure 11 
(left panel) shows the dispersion relation in such a plasma and the resonance conditions for 
electrons and protons with (3 = 0.5 and // = 0.12. 

With this additional branch, particles with /j, ^ can interact at least with three waves 
with one from each of the three cyclotron branches (particles with /i — always interact 
with two waves from one of the wave branches). Some particles can resonate with five waves 
three of which would be from one of the three cyclotron branches (e.g. the protons in the 
left panel of Figure 11 have three resonances with the EC branch with two of them shown 
in the lower panel of the figure. The third one is at high k and w — > Q c ). In a strongly 
magnetized plasma the two additional resonances can also come from the electromagnetic 
branches. The general results are quite similar to those in pure hydrogen plasmas; there 
are critical angles and critical energies, which now separate the /i — (5 and (3 — a spaces 
into regions with three and five resonant interactions. Given spectra for each of the wave 
branches one can proceed to calculate the F-P coefficients, which have sharp jumps across 
the critical angles or energies, and the times r ac and r sc . In the next section we discuss the 
turbulence spectrum we shall use for this purpose. 



-33- 



c; 

\ 

3 



a 

\ 

3 



1 1 1 1 

^- ^ a = 1.0 


Y„i = 0-08^^ ■ 


- 




M = 0.12 

1 


(3 = 0.5 

1 



kc/n e 




n 1 1 1 1 1 1 1 1 1 1 — 

q = 4.0-M.7, a = 0.98 




Electron 



_L 



_L 



_L 



2 4 6 

log 10 (E/keV) 



Fig. 11. — Left panel: Same as Figure 1 but in an electron-proton- 4 He plasma with a = 1.0 and 
Yfie = 0.08. The resonant interactions (circles) are for electrons (upper line) and protons (lower 
line) with (5 = 0.5 and fi = 0.12. Here we note that there are two ion cyclotron branches (both 
indicated by the dotted curves), the upper one approaches to the He cyclotron waves at high k 
(HeC) while the lower one approaches to the modified proton cyclotron waves (PC). The crosses 
indicate the expected breaks in the turbulence spectrum as discussed in the text. Right panel: 
The pitch angle averaged acceleration (solid curves) and scattering (dashed curves) times in units 
of r p for electrons (thick curves) and protons (thin curves). The model parameters are indicated in 
the figure. The iiiiniiiiLiiii of the electron acceleration time corresponds to the spectral break & max 
of the turbulence in the EC branch. See text for details. 
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4.2. Spectrum of Turbulence 

In the discussion of the previous section we assumed a power law wave spectrum £ oc k~ q , 
for k > /c mm , with fc m i n <C k\ = ad 1 ^ 2 so that there will be waves with the Alfvenic dispersion 
relation for interactions with high energy electrons and protons. We now introduce an 
additional cutoff at high wavenumbers presumably caused by thermal damping. Thus we 
will have a broken power-law turbulence spectrum with three indexes q, qi, and % and two 
critical wavenumbers k m \ n and A; max - 

{{k/k m \ n ) qi ) for k < /Cmin! 

(k/k miQ )- q , for k min < k < k max ; (29) 

(^max/^min) ? (^/^max) 9h ) f° r k > & m ax; 

where q\ > (we choose q\ — 2 because its value is almost irrelevant), g = 1.7 is the 
Kolmogorov value, and q^ = 4, a typical value of the spectral index for waves subject to strong 
damping (Vestuto, Ostriker & Stone 2003). A self-consistent treatment of wave-particle 
interactions is required for an exact description of these spectral breaks. This is beyond the 
scope of the current investigation. In stead, we make some reasonable assumptions on these 
breaks: 

The low k or large scale cutoff &; mm > c/(f2 e L max ) where L max is the largest scale of 
the turbulence, which must be less than the size L of the region, and is most likely much 
less than it. To accelerate particles to high energies, we also need k min <C A; a — ex/ A3. For 
most waves we choose k m i n as before, i.e. by the value of the highest energy we want the 
particles to achieve. However, such a choice for the PC branch results in a sharp feature in 
the spectrum of the accelerated protons. As can be seen in the lower left panel of Figure 11, 
the PC branch, unlike the EC, PC or HeC branches, crosses the frequency axis (k = 0) at 
a finite frequency wpc 7^ 0. Such waves can resonantly scatter protons with a Lorentz factor 
of ~ 2/(1 + 2Yhc) ~ 1-7 or an energy of a few hundreds of MeV. If the spectrum of PC 
branch wave extends to a very low wavenumber, one would get very efficient acceleration at 
such energies and a sharp spectral feature. We assume that such a feature is not present 
(although there is no definite observation to rule it out) and cutoff the spectrum of the PC 
waves at a higher k niin ~ A:a/5 or a scale of L max ~ 4007rc/o;fi c . 

The high k or small scale cutoff is determined via the damping of the waves by low energy 
particles. For example, the cyclotron waves with high wavenumbers are subject to thermal 
damping in plasmas with a finite temperature (Schlickeiser & Achatz 1993; Steinacker et al. 
1997; Pryadko & Petrosian 1998). One can introduce an imaginary 00 to include this effect. 
The real part of uo is not very sensitive to the plasma temperature (Miller & Steinacker 
1992). So the cold plasma dispersion relation still gives a good description of the resonant 
interactions between waves and particles. However, above a wavenumber where the thermal 
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damping time is comparable with the time scale of the wave cascade, the thermal damping 
steepens the spectrum of the turbulence. In the absence of a full treatment of these processes, 
we shall assume that the cyclotron waves have steeper spectra than the Whistler and Alfven 
waves and set A; max — — a for the EC branch. For the PC and HeC branches we set 
fc max = k\ = a5 l l 2 . These spectral breaks are indicated by the crosses in the lower left panel 
of Figure 11. 

Recent studies of the transport of high energy particles in the solar wind indicate that 
the resonance broadening due to the dissipation of turbulent plasma waves plays an essential 
role in explaining the observed long mean free paths of the particles (Bieber et al. 1994; Droge 
2003). Such broadening is also expected in our case. It will modify only equation (6) if the 
broadening width is less than the separation between the resonances. In the opposite (and 
unlikely) case the idea of resonant interaction is invalid. Because the resonance broadening 
mostly affects the scatterings of particles at relatively low energies, where Coulomb collisions 
are important under solar flare conditions, it is less important in understanding the particle 
acceleration processes studied here. 

All the formulae developed in the previous sections are still valid except that now 

(9 " , (30) 

and the Alfven velocity is given by Pa = S 1 ^ 2 /[a(l + 2YH e ) 1 / 2 ]. Because q\ > and > q 
(and the PC branch contains much less energy than the other branches) the total turbulence 
energy density £ tot ~ £ Q . 

The right panel of Figure 11 gives the electron (thick curves) and proton (thin curves) 
acceleration (solid curves) and scattering (dashed curves) times for a model with a — 1. 
We see that in the high energy range where particles are mostly interacting with the Alfven 
waves, the times are almost the same as those in a pure hydrogen plasma. At low energies, the 
times rise sharply with the decrease of energy due to the thermal damping of the waves with 
high wavenumbers. As we will see in the following discussion, the thermal damping makes 
the particle acceleration times match their escape times, giving near power law accelerated 
particle distributions. More importantly, in the intermediate energy range, an acceleration 
barrier in the proton acceleration time still exists even though it is not as prominent as 
it is in a pure hydrogen plasma. This is mainly due to interactions with the HeC branch 
which makes the acceleration of low energy protons more efficient. The sharp decrease of 
the acceleration time with energy near the critical energy is still due to interaction with 
the Whistler waves. Comparing with Figure 8, we note that the electron acceleration and 
scattering times are also affected as evident by the wiggles seen at a few tens of MeV. These 
wiggles are due to interactions with the HeC and PC branches. 
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We emphasize here that q = 1.7 corresponds to the Kolmogorov spectrum. Our only 
assumption which may be ad hoc is the large scale size cutoff for the PC branch waves. 
This assumption is not driven by observations but primarily introduced to obtain a smooth 
proton spectrum. We shall explore consequences of the assumption in the future. In the 
following discussion, we will fix these parameters at the specified values and investigate how 
the particle acceleration processes are affected by the strength of the turbulence £ tot , the 
size of acceleration site L, the plasma parameter a and the temperature of the injected 
particles. We will show that this model gives much more reasonable explanations to solar 
flare observations than the previous one. 



4.3. Relative Acceleration of Electrons and Protons 

We now present some results on the relative numbers of accelerated electrons and protons 
at the acceleration site (LT) and escaping to the FPs. Here we explore its dependence on the 
model parameters. The normalization is as before (see Figure 9 and discussion in eq. [21]), 
namely we assume that the escape fluxes for the total numbers of electrons and protons are 
equal. 

Figure 12 gives our fiducial model for solar flares where the energy content in the 
accelerated electrons and protons in the relevant observational energy bands (indicated by 
the shaded regions) is comparable. The time scales are given in the left panel and the 
corresponding particle distributions are shown in the right panel. The line types are the 
same as those in Figure 9. The temperatures of the injected electrons and protons are 
T = 1.5 keV. The other model parameters are: the size of the LT source L = 5 x 10 s 
cm, n e = 1.5 x 10 10 cm~ 3 , B = 400 Gauss (a = 0.98) and r p = 1/70 sec. These imply 
87r£ to t/5o - 8 - 8 x lO -9 ^ 7 , which is < 10~ 3 for k min > 2nc/Ltt e = 5.2 x HT 8 . Compared 
with models in pure hydrogen plasmas, the turbulence required to accelerate the particles 
is much weaker in the current model. This is mainly because of the adopted Kolmogorov 
(instead of q = 2 or 3) turbulence spectrum at long wavelengths. 




Fig. 12. — Same as Figure 9 but for different model parameters (listed) and in a hydrogen and 
helium plasma. The energy content in high energy electrons and protons is comparable. The 
hatched regions in the right panel correspond to the energy bands related to observations of the 
hard X-ray and gamma-ray emissions during the impulsive phase of solar flares. The signs indicate 
the spectral breaks of the accelerated particle distributions. 
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log 10 (E/keV) log 10 (E/keV) 



Fig. 13. — Same as the right panel of Figure 12 but for two different values of r p , one larger (left 
panel) and one smaller (right panel). 
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In Figure 13, we show the effects of the strength of turbulence on the particle accel- 
eration. Because the acceleration time is proportional to r p but the escape time is inverse 
proportional to it, a small change in t p can make the acceleration and escape times off bal- 
ance very quickly. The spectra of the accelerated particles then change dramatically with r p . 
However, because this is true for both electrons and protons, changes in the energy partition 
between electrons and protons are much smaller. 

Another parameter that affects the spectra of the accelerated particles is the size L of 
the acceleration region. An increase of L results in an increase of the escape time and harder 
spectra of the accelerated particles. However, again because the escape times of electrons 
and protons increase by the same factor, the energy partition between them is not changed 
significantly. For example, for a model with L = 10 9 cm and r^ 1 = 50 s" 1 , we find that the 
accelerated particle spectra which are harder than those in the left panel of Figure 13 and 
more like those in the right panel of Figure 12. 

Next we examine the effects of the plasma parameter a oc n 1 / 2 /B. Figure 14 shows 
how the relative acceleration of electrons and protons changes with the change in the value 
of a. It turns out that for the range of parameters used in the current study, it does not 
matter whether a is changed by changing the value of the density n or the magnetic field 
B. The difference between these two possibilities will appear as relatively small changes 
in the spectra at the low and high energy ends where Coulomb collisions and synchrotron 
losses become important, respectively. To make the spectral shapes compatible with solar 
flare observations, the value of r p also needs adjustment. But as we showed above, r p affects 
primarily the spectral hardness but not the relative acceleration of electrons and protons. 
Thus the most relevant cause of the changes in the relative acceleration of protons and 
electrons is the variation of a; proton (and consequently other ions) acceleration is more 
efficient in high density and/or low magnetic field plasmas. Given that the acceleration of 
electrons and protons are dominated by different wave branches, it is not surprising that 
their relative acceleration depends on the plasma parameter a. 
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Fig. 14. — Same as Figure 12 but for models with different plasma parameters a, which are 
indicated in the figures. r p is chosen such that the accelerated particle distributions are similar 
to that required to explain the impulsive hard X-ray emission from solar flares. The other model 
parameters are the same as those shown in Figure 12. We see that the electron acceleration 
is favorable in strongly magnetized plasmas while the proton acceleration dominates in weakly 
magnetized plasmas. 
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Fig. 15. — Same as Figure 12 but for models with different temperatures for the injection plasma. 
The temperatures are indicated in the figures. Other model parameters remain the same as the 
fiducial model. The distributions of high energy particles are not affected by the injection process. 
However, at low energies where Coulomb collisions dominate, much less particles are accelerated 
for an injection plasma with a lower temperature. 
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Finally, we consider the effects of the temperature of the background or injected par- 
ticles. In the models discussed above, we use a high value of temperature of a few keV, 
which requires a pre-heating of the flaring plasma to a temperature above the quiet coro- 
nal value. GOES and RHESSI observations do suggest such a characteristic energy for the 
particles before the impulsive phase of X-ray flares. For example, RHESSFs high resolution 
spectra indicate that the electrons in the soft X-ray emitting plasma always deviate from 
an isothermal distribution, implying a significant pre-heating of the flare plasma (Holman 
et al. 2003). The effects of the injection temperature are demonstrated in Figure 15, which 
shows particle spectra for models with temperatures different than that in the fiducial model 
(Figure 12). All other model parameters remain the same. We see that the shapes of the 
spectra in the high energy range do not change significantly. However, with the increase of 
the temperature, more particles reach the energy range where the acceleration rate is larger 
than the Coulomb collisional loss rate and are eventually accelerated to higher energies. At 
lower temperature, the quasi-thermal part of the spectra (similar to that of the injected 
particles) is more prominent, while at higher T the spectra of the accelerated particles are 
dominated by the nonthermal tails. 



5. SUMMARY AND DISCUSSION 

The primary aim of this work is the determination of the relative acceleration of electrons 
and protons by waves in a stochastic acceleration (SA) model. In this paper, we present the 
results of the investigation of the resonant interaction of the particles with a broad spectrum 
of waves propagating parallel to the large scale magnetic field. We calculate the acceleration 
and transport coefficients and determine the resulting spectra for both particles for physical 
conditions appropriate for solar flares. We show that the injection of such a turbulence in a 
magnetized hot plasma can accelerate both electrons and protons of the thermal background 
plasma to high energies in the acceleration site. Some of the accelerated particles escape 
the site and reach the FPs. The parameters that govern these processes are the density, 
temperature, magnetic field, size of the acceleration region and the intensity and spectrum 
of the turbulence. 

We first describe two general features of our results. 

A. The first has to do with the general characteristics of the accelerated particle spectra. 
The outcome of the SA of a background thermal plasma is the presence of two distinct com- 
ponents. The first is a quasi-thermal component at low energies where Coulomb collisions 
play important roles. This can be considered as a simple heating process of the background 
plasma. The second is a nonthermal tail with a somewhat complex spectral shape. Techni- 



-43- 



cally, one can separate the two components at the energy where the Coulomb collisional loss 
rate -Ecoui is equal to the direct acceleration rate A. We can then calculate the fractions of 
the turbulence energy that go to "heat" and to "acceleration" . This explains the observation 
of both thermal and nonthermal emissions during the impulsive phase of solar flares. 

The spectra of particles reaching the FPs are in general harder than the corresponding 
LT particle spectra because high energy particles in the nonthermal tails escape more readily. 
The relative size of the two components in the acceleration site (LT) vs the FPs and for 
electrons vs protons depends sensitively on the model parameters, which can explain the large 
variation of the observed nonthermal emission among solar flares and other astrophysical 
sources. 

B. The second feature has to do with the relative acceleration of electrons and protons. 
To our knowledge, a new result of our investigation is that there is a significant difference in 
the acceleration of protons and electrons. While the transport and acceleration coefficients 
for electrons are smooth functions of energy, this is not true for protons. There appears to 
be a barrier (lower rate) of acceleration for intermediate energy protons. This can have a 
dramatic effect on the relative production rate and spectra of the accelerated protons and 
electrons. 

To demonstrate some of this and other more subtle effects, we first investigated the 
acceleration of electrons and protons in pure hydrogen plasmas by turbulence with a simple 
power-law spectrum. We find that this simple model does not agree with some qualitative 
aspects of the observed accelerated particle distributions in solar flares. The barrier for pro- 
tons is too strong for reasonable physical conditions. We then explore more realistic models, 
where we include the effects of the background 4 He particles and the thermal damping of the 
waves. These more realistic models are in better concordance with solar flare observations. 

Specifically we find the following: 

1. In general, electrons are preferentially accelerated in more strongly magnetized plasmas 
(small a oc n 1//2 / B) while the proton acceleration is efficient in more weakly magnetized 
plasmas. The ratio of the energy that goes into the accelerated electrons to that 
into protons is very sensitive to a, which can explain the wide range of the observed 
energy partition between these particles. The proton acceleration will be more efficient 
in larger loops where the magnetic field is presumably weaker and during the later 
phase of flares when the corona loops have been filled by plasmas evaporated from the 
chromosphere, giving a higher gas density. This can explain the offset of the centroid 
of the gamma-ray line emission (due to accelerated protons and ions) from that of the 
hard X-rays indicated by a recent RHESSI observation (Hurford et al. 2003). It can 



-44- 



also account for the observed delay of the nuclear line emission relative to the hard 
X-ray emission (Chupp 1990). 

2. The acceleration rates and spectra of both electrons and protons are very sensitive to 
the intensity of the turbulence and the size of the acceleration site. Models with more 
intense turbulence and/or larger acceleration region give rise to harder spectra. This 
result can explain the observed soft X-ray emission in advance of the impulsive phase 
hard X-ray and gamma-ray emissions (the so-called pre-heating) and the slower than 
expected decline of the temperature of the LT plasma in the gradual phase. When the 
turbulence is weak, as will be the case at the beginning and end of a flare, almost all 
the dissipated turbulence energy goes into the quasi-thermal component and there is 
no significant hard component, producing soft X-ray emission without obvious hard X- 
ray or gamma-ray emission. When the strength of the turbulence exceeds a threshold, 
nonthermal tails and high energy radiations ensue. On the other hand, for a turbulence 
energy much above this threshold, one would expect harder spectra than observed in 
solar flares. This may indicate that the sudden presence of a large amount of high 
energy particles also introduces significant dissipation of the turbulence over a broad 
frequency range such that the strength of the turbulence is limited to a level close to 
the threshold. Consequently, we do not see flares with very flat X-ray spectra. To 
address these processes in detail, one needs to treat the wave generation, cascade and 
damping by both low and high energy particles properly. Such an investigation is 
clearly warranted now but is beyond the scope of this paper. 

3. In general the spectra of both electrons and protons at the acceleration site (LT) are 
softer (stronger quasi-thermal component and steeper nonthermal spectrum) than the 
equivalent thick target spectra at the FPs. This is in excellent agreement with the 
results from the YOHKOH (Petrosian et al. 2002) and with the more convincing evi- 
dence from RHESSI observations (Jiang et al. 2003). The most important parameter 
here is the energy dependence of the escape time (see equation [16]) which depends on 
the pitch angle diffusion coefficient and the size of the acceleration region. Unlike the 
acceleration time (see item B above), the scattering and escape times for protons and 
electrons have similar general behaviors. Consequently, the difference between the LT 
and FP spectra is similar for both electrons and protons. 

4. For injected plasmas with high temperatures, most of the particles can be accelerated 
to very high energies and the steady state particle distribution at low energies can be 
quite different from a thermal distribution. The presence of a quasi-thermal component 
is typical for low temperature plasmas. 

5. There are high energy cutoffs (at around 1 MeV for electrons and 10 MeV for protons) 
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in the accelerated particle spectra. Both cutoffs are due to the quick escape and the 
relatively inefficient acceleration of higher energy particles. The location of these cut- 
offs are directly related to the higher wavenumber spectral breaks in the turbulence 
spectrum. In plasmas with stronger thermal damping, the acceleration of high energy 
particles becomes relatively more efficient than that of low energy particles. Conse- 
quently, the cutoffs shift toward higher energies. We would then expect a positive 
correlation between the cutoff energies and the heating rate of the background plasma. 
Observations over a broad energy range will be able to test this prediction. 

Finally, we summarize several improvements that are required for direct comparisons 
with observations: 

1. The results presented here show that wave-particle interactions play crucial roles in 
solar flares, especially during the impulsive phase. A self-consistent treatment of this 
problem requires the solution of the coupled kinetic equations for both particles and 
waves. Previous studies on this aspect focused on the Alfven waves alone and ignored 
the energy dependent escaping processes (Miller & Roberts 1995; Miller et al. 1996). 
Incorporation of the current investigation will make the models more realistic. 

Waves propagating obliquely with respect to the large scale magnetic field will intro- 
duce new features to the wave-particle interaction (Pryadko & Petrosian 1999). Earlier 
studies have shown that the fast-mode waves are very efficient in heating or accelerat- 
ing super- Alfvenic particles via Landau damping or transit-time damping (Miller et al. 
1996; Quataert 1998; Schlickeiser & Miller 1998). These waves are expected to enhance 
the acceleration of high energy electrons and protons and sub-Alfvenic particles may 
also be accelerated when one adopts the exact dispersion relation for the waves. Results 
similar to what we present here are expected. Moreover, if the turbulence is dominated 
by the lower hybrid waves, electrons will not be scattered efficiently so that the electron 
distribution is not isotropic (Luo, Wei, & Feng 2003). The acceleration barrier may 
not exist if this is also true for protons. A comprehensive study including these waves 
is needed to address the heating and acceleration processes more completely. 

2. A time dependent model is needed to address the temporal characteristics of solar 
flares and the injection processes. Here we assume that the system is in a steady 
state and the injection fluxes of protons and electrons are equal. This may be the 
case if the plasmas are brought into the acceleration site by the reconnecting magnetic 
fields. However, e.g. if electrons have a shorter escape time than protons, there could 
be a net charge flux from the acceleration site, which would induce reverse currents 
consisting mainly of electrons so that the injection fluxes of electrons and protons into 
the acceleration site will be different. 
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3. The application of the formalism developed here to the acceleration of other ions is 
straightforward. We are in the process of evaluating the relative acceleration of different 
ion species and isotopes and the results are promising and will be published in future 
papers. It is also straightforward to apply the formalism to accretion systems of black 
holes and neutron stars. Besides the magnetic reconnection, turbulent plasma waves 
can also be produced by the magneto-rotational instability in accretion disks (Balbus 
& Hawley 1991). 



The work is supported by NASA grants NAG5-12111, NAG5 11918-1, and NSF grant 
ATM-0312344. 



A. Dispersion Relations for the EC and PC Branches 

Waves in the EM' branch interact resonantly with protons and those in the EM branch 
interact with electrons only for low values of a. These interactions mostly affect the acceler- 
ation of very low energy particles (see discussion in § 3.2). Waves in the EC and PC branches 
are the dominant modes for the acceleration of protons and electrons for intermediate and 
high values of a and energies and therefore play key roles in determining the relative accel- 
eration of the two species. In what follows we give some approximate analytic descriptions 
of these modes which are considerably simpler than equation (22). 

At frequencies w < 5, or \k\ <C &a = ad 1 ^ 2 , both branches reduce to the Alfven 
waves with the dispersion relation uj = \k\f3\ — 8{\k\/k\) for the EC branch. The middle 
portion of the EC branch, k& <C \k\ <C k w = a, corresponds to the Whistler waves with 
the dispersion relation uo ~ k 2 /a 2 = (/c//c w ) 2 . At still higher wavenumbers, (\k\ 3> few), 
uj — > 1 and the dispersion relation can be approximated as uj ~ 1 — a 2 /k 2 . The transition 
between the Whistler and this portion occurs at few = oc. This suggests that the dispersion 
relation for the Whistler and electron-cyclotron portions can be approximately described by 
uj = k 2 /(k 2 + a 2 ). For the EC branch we can use the simple approximation 

f) l l 2 X + x 2 

uj — — — ; with x = \k\/k w = \k\/a , (Al) 

1 + x z 

which agrees with the exact expression within 40% for a > 0.6. For highly magnetized 
plasmas with a < 1 the Whistler branch disappears. One then has uj — \k\ for \k\ < 1 and 
uj = 1 for \k\ > 1 for the EC branch (and the reverse is true for the EM branch). 

Similarly for the PC branch one gets the Alfven waves with uj = — \k\0A at low wavenum- 
bers (\k\ <IC &a). The proton-cyclotron waves whose dispersion relation doesn't have a simple 
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form can be roughly approximated as Upc ~ —5(1 — 1/(2 + (k/kx) 2 ) for high k. We can 
combine these two forms into one simple expression 

" = -a /tLa ; with = I*IAa = |*|/(«* 1/2 ), (A2) 

which agrees with the exact expression to within 10% for a > 0.6. For very highly magnetized 
plasmas, a <C <5 1//2 , one again has uo = —\k\ for \k\ < 5 and uo = — 5 for \k\ > 5 for the PC 
branch (and the reverse for the EM' branch). 



af3x = — — — - — , x = k/kw = k/a. (Bl) 



B. Critical Energies and Angles for Resonance with the EC and PC Branches 

With the approximate analytical expressions (Al) and (A2) for the dispersion relation, 
one can derive the critical velocity (24) and critical angle (25) for resonant interactions of 
low energy particles with the EC and PC branches. The critical velocity is for protons with 
/i — 1 interacting with the EC branch. From the resonance condition (4) and the dispersion 
relation (Al) for the EC branch, we have 

5 _ 5 1/2 \x\ +x 2 
7 1 + x 2 

This equation has three roots with two of them being equal and <C 1 at the critical velocities 
(see Figure 1, and equation [24]). We can therefore ignore the x 2 term in the denominator 
of the right hand side of equation (Bl). One can show that 

p CI = 3V5/a ~ 0.07/a , (B2) 

which agrees with the numerical result within 15% [eq.(24)]. In general, we have 

p CT = 3V5/fia ~ 0.07/^a . (B3) 

For electron resonances with the EC branch, we have 

1 _ 5_V 2 \x\+ x 2 
7 1 + x 

The equation has three roots with two of them being equal when 



aPfix + - = — 17~^2 — • ( B4 ) 



(a/37/icr) 2 = ; 8(7 • (B5) 

8 + 20 7 - 7 2 + y/{8 + 20 7 - 7 2 ) 2 + 64( 7 - l) 3 

We then have 

H a ~P 2 /V54a, for (3 < 1 , (B6) 
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which agrees with equation (25) within 5%. 

For proton resonance with the PC branch, we have 

£ I I i 2 

aP/iV5y -- = -6 - |y \ + V 5 , y = k/V5a. (B7) 

For low energy protons at critical pitch angles, the two equal roots of the equation are much 
larger than unity. We can therefore approximate the right hand side of this equation as 
—5(1 — l/y 2 )- Then we have 

4 5 _ ( 7 
27 ~ Vt 3 ! 

which becomes 

Atcr = y/5p 2 /y/5Aa for /3 < 1 . (B9) 
This agrees with equation (25) within 5%. 




C. Approximate Analytic Expressions for the Acceleration and Scattering 
Times in the Relativistic and Low Energy Limits 

It is useful to have some approximate analytical expressions for the acceleration and 
scattering times under certain limits. In the relativistic region where 7 ^> 1, the results for 
electrons are relatively simple and have been studied under different context (Schlickeiser 
1989; Pryadko & Petrosian 1997). Here we discuss the results for ions. For relativistic 
particles, and in general for weakly magnetized plasmas (a > 5 1 / 2 ), R\ ^> As a result, 
the acceleration time defined by equation (12) can be approximated as: 



r. 



= (D PP / P 2 ) , (CI) 



where "<>" denotes average over pitch angle. Relativistic particles with 7 \u>i\/5 resonate 
with the Alfven waves with u = ±j3\k. From the resonance condition (4), we have the 
wavenumbers of the resonant waves: k± = Ui/[y((3[i =F /3a)] — Uififx/j. Then we have 

D PP T P i \ = 7 g 2 Pl f 1 A n h_ n 2\\na.,_a\\q-l,(\a.,,a\\q-l] 



l^-^-^^M^' (C2) 



V 2 I H\ q f3 2 

„ f^K C ..2,,, , „, , i-,"- 2 i"- : M 

\Ui\"P 2 



which is consistent with the numerical results within a factor of two. 
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Similarly, one can estimate the scattering time in the relativistic limit: 



Tpi \ (P/i - 



(C3) 



2-q A-q 

which is in agreement with the numerical results within 50%. Note that the integral over \i 
has been taken from f3\/ (3 to 1 here. Like the acceleration time, protons and electrons with 
the same energy have the same scattering time. This is consistent with the result of Pryadko 
& Petrosian (1997). We note that in the relativistic region, the acceleration and scattering 
times are identical for all charged particles except for the Ui term. 

In the nonrelativistic region where 7 ~ 1, from the resonance condition uo = —Ui + [3fik, 
and the dispersion relation (22), one can show that 

k = -ua{(3fik)- 1/2 . (C4) 

So k = (ujfa 2 ^' 1 ^ 1 ) 1 / 3 . Using equation (5), we have 



D r 

^pp 1 pi 

p2 




2 2\ (1-9) /3 V 



-2( g +2)| ,2(1 



-9)^-7] 



71 l/3 



(C5) 



(g + 2)(g + 8) 

—2(3fi under the resonance condition. In a previous study (Pryadko & 
Petrosian 1997), the minus sign was missed, which causes their acceleration time three times 
shorter than ours. Equation (5) then gives 



We note that j3„ 



1 + 



(x(\ui\af3fx) 



1/3 



Pa 



u?a 2 



uJi a 



-2q/3 



1 



-9/3 > 



2bJ 1 /3 



g(? + 6) (g + 4)(g + 10)(a/5) 2 / 3 



a;, 



1 2/3 



(q + 8)(q+U)(af3)^\ ' 
which agrees with the numerical results within 20%. 



(C6) 



D. Approximate Analytic Expression for the Acceleration Barrier 



To estimate the acceleration time at the barrier in pure hydrogen plasmas, we notice that 
resonant interactions of protons with nearly 90° pitch angle (/i < // cr ) with proton-cyclotron 
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waves moving in both directions have significant contributions to the proton acceleration 
below the critical energy. Because D^(Ri — R 2 ) is a smooth function near /i = (Figure 5) 
and /x cr <C 1 (eq. [25]), we have 



fac(p) 



(Dl) 



Beyond the critical energy, some protons start to resonate with Whistler waves (see Figure 
2) and the proton acceleration rate increases sharply with energy. 



From equation (5), we have 



D pp ( f i = 0,p)_ 2(3 2 ph \k\ 



-q 



wrw g i ' 



(D2) 



where k can be obtained from the dispersion relation (22) and u = — 5/ 7 is given by the 
resonance condition (4). Then we have 



1 + 



gV(i+j) 

6(6 + i)(t-1) 



2/i cr 

0(9-l)/2 T X(q-5)/2 q+1 
A*cr 

= 7a' ?+2 (25)^- 1 )/ 2 /? 1 -^- 5 / 2 r p 
7a q+2 6 "' 2 / 



(3+?)/2 



1 | a 2 (l + £){ 7 - 0.5(1 -£)} 7 



(5 + 7 ) 2 ( 7 -l) 2 <5 



for /3 < 1 



'p ■ 



(D3) 



Combining this with equations (B2) and (C2), one can estimate the height of the acceleration 
barrier in logarithmic scale at the critical velocity 



51og(r ac ) ^ log {3.1 x 2^- 1 )/ 2 a« +2 r 1 / 2 -V[g(g + 2)]}. 
Similarly, one can estimate the acceleration time for low energy electrons: 



(D4) 



r P 7 2 -^ 2 

2/icr 

7a 9+2 



1 + 



E 



m c c z 



a 2 ^ 2 (l + 5) 
( 7 -l)(l + 7 (J) 

(l-5)/2 



(3+g)/2 



1 + 



a 2 (l + <5){ 7 5 + 0.5(1 -5)} 7 3 



( 7 _1)2( 1+7j5 )2 



for /3 < 1 



These expressions are consistent with the numerical results within a factor of two. The 
discrepancy is large for turbulence with a flat spectrum. This is mainly due to contributions 
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from electron-cyclotron and proton-cyclotron waves to the acceleration of particles with two 
resonances. When the turbulence spectrum becomes flatter, their contributions to the pitch 
angle averaged acceleration time becomes more important. However, as we discussed in 
§ 4, this effect is not important in real astrophysical situation where the cyclotron waves 
are damped. So the analytical expressions give a good estimate of acceleration time in the 
intermediate energy range. 

REFERENCES 

Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 
Barbosa, D. D. 1979, ApJ, 233, 383 

Bech, F. W., Steinacker, J., & Schlickeiser, R. 1990, Solar Phys. 129, 195 

Bieber, J. W. et al. 1994, ApJ, 420, 294 

Chupp, E. L. 1990, Science, 250, 229 

Drdge, W. 2003, ApJ, 589, 1027 

Dung, R., & Petrosian, V. 1994, ApJ, 421, 555 

Ginzburg, V. L., & Syrovatskii, S. I. 1964, The Origin of Cosmic Rays, (Pergamon Press, 
New York) 

Hamilton, R. J., k Petrosian, V. 1992, ApJ, 398, 350 

Holman, G., Sui, L., Schwartz, R., & Emslie, G. 2003, ApJ, 595, L97 

Hua, X. M., & Lingenfelter, R. E. 1987a, ApJ, 323, 779 

Hua, X. M., & Lingenfelter, R. E. 1987b, Sol.Phy., 107, 351 

Hua, X. M., Ramaty, R., & Lingenfelter, R. E. 1989, ApJ, 341, 516 

Hurford, G. J., Schwartz, R. A., Krucker, S., Lin, R. P., Smith, D. M., & Vilmer, N, 2003, 
ApJ, 595, L77 

Jiang, Y. W., Liu, W., Petrosian, V., & McTiernan, J. 2003, AAS/SPD, #18.03 
Kennel, C. F., & Engelmann, F. 1966, Phys. Fluids, 9, 2377 
Kirk, J. G., Schneider, P., & Schlickeiser, R. 1988, ApJ, 328, 269 



-52 - 



Kozlovsky, B., Murphy, R. J., k Ramaty, R. 2002, ApJS, 141, 523 
Leach, J. 1984, Ph.D. thesis, Stanford Univ. 

Liu, W., Jiang, Y. W., Petrosian, V., k Metcalf, T. R. 2003, AAS/SPD, #34, #18.02 
Luo, Q. Y., Wei, F. S., k Feng, X. S. 2003, ApJ, 584, 497 
Masuda, S. 1994, Ph.D. thesis, Univ. Tokyo 

Masuda, S., Kosugi, T., Hara, H., Tsuneta, S., k Ogawara, Y. 1994, Nature, 371, 495 

Mazur, J. E. Mason, G. M., Klecker, B., k McGuire, R. E. 1992, ApJ, 401, 398 

Miller, J. A. 2003, Multi- Wavelength Observations of Coronal Structure and Dynamics eds. 
Petrus C.H. et al. (COSPAR Colloquia Series Vol. 13), 387 

Miller, J. A., Cargill, P. J., Emslie, G., Holman, G. D., Dennis, B. R., LaRosa, T. N., 
Winglee, R. M., Benka, S. G., k Tsuneta, S. 1997, J.Geophys.Res. 102, 14631 

Miller, J. A., LaRosa, T. N., & Moore, R. L. 1996, ApJ, 461, 445 

Miller, J. A., & Ramaty, R. 1987, Sol.Phy. 113, 195 

Miller, J. A., & Roberts, D. A. 1995, ApJ, 452, 912 

Miller, J. A., & Steinacker, J. 1992, ApJ, 399, 284 

Park, B. T., k Petrosian, V. 1995, ApJ, 446, 699 

Park, B. T., k Petrosian, V. 1996, ApJS, 103, 255 

Park, B. T., Petrosian, V., k Schwartz, R. A. 1997, ApJ, 489, 358 

Petrosian, V., k Donaghy, T. Q. 1999, ApJ, 527, 945 

Petrosian, V., Donaghy, T. Q., & McTiernan, J. M. 2002, ApJ, 569, 459 

Petrosian, V., McTiernan, J. M., & Marschhauser, H. 1994, ApJ, 434, 747 

Post, E. F. 1956, Reviews of Modern Physics, 28, 338 

Pryadko, J. M., & Petrosian, V. 1997, ApJ, 482, 774 

Pryadko, J. M., k Petrosian, V. 1998, ApJ, 495, 377 

Pryadko, J. M., k Petrosian, V. 1999, ApJ, 515, 873 



-53- 

Quataert, E. 1998, ApJ, 500, 978 

Ramaty, R. 1979, in Particle Acceleration Mechanisms in Astrophysics, eds. J. Arons, C. 
Max, k C. McKee, (New York: AIP) 

Ramaty, R., Kozlovsky, B., k Lingenfelter, R. E. 1975, SSRv, 18, 341 

Ramaty, R., Kozlovsky, B., k Lingenfelter, R. E. 1979, ApJS, 40, 487 

Reames, D. V., Meyer, J. P., k von Rosenvinge, T. T. 1994, ApJS, 90, 649 

Rieger, E., Gan, W. Q., k Marschhauser, H. 1998, Sol.Phy. 183, 123 

Schlickeiser, R. 1989, ApJ, 336, 243 

Schlickeiser, R. 2002, Cosmic Ray Astrophysics (Berlin: Springer) 
Schlickeiser, R., k Achatz, U. 1993, J. Plasma Phys., 49(1), 63 
Schlickeiser, R., Miller, J. A. 1998, ApJ, 492, 352 

Spitzer, L. 1956, Physics of Fully Ionized Gases (Interscience Publishers Inc. New York) 

Steinacker, J., Meyer, J. P., Steinacker, A., k Reames, D. V. 1997, ApJ, 476, 403 

Steinacker, J., k Miller, J. A. 1992, ApJ, 393, 764 

Sturrock, P. 1994, Plasma Physics (Cambridge University Press) 

Vestuto, J., Ostriker, E., k Stone, J. M. 2003, ApJ, 590, 858 



This preprint was prepared with the AAS IATgX macros v5.2. 



